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0. Introduction 

The study of the spectrum of coupled random matrices has received rather 
little attention. To the best of our knowledge, coupled random matrices have 
been studied, to some extent, by Mehta in [16], [17], [11]. In this work, we ex- 
plain how the integrable technology can be brought to bear to gain insight into 
the nature of the distribution of the spectrum of coupled Hermitean random 
matrices and the equations the associated probabilities satisfy. In particular, 
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the two-Toda lattice, its algebra of symmetries and its vertex operators will 
play a prominent role in this interaction. Namely, the method is to introduce 
time parameters, in an artificial way, and to dress up a certain matrix inte- 
gral with a vertex integral operator, for which we find Virasoro-like differential 
equations; for a state of the art survey of these methods, see [20]. These meth- 
ods lead to very simple nonlinear third-order partial differential equations for 
the joint statistics of the spectra of two coupled Gaussian random matrices. 

Bi-orthogonal polynomials and the two-Toda lattice. Given a weight, de- 
pending on parameters t = (ti, £2, • • •) and s = (s±, S2, ■ ■ •) £ C°° and a coupling 
constant c, 

,.,....„, .., „ . 1 iX +cx y-^i s iv ; 

consider the associated (monic) bi-orthogonal polynomials pp (x) and p\^ (y) 
of degree i defined by 

(Pk^pf^) = h kSk,e, where (f,g) = / / dxdy f(x)g(y)p(x,y). 

J JR 2 

In [2], it was shown that the associated pair of semi-infinite vectors, depending 
on (t,s), 

$ x ( z ) : = eEr^V 1 )^) and ^*(z) := e'^T SkZ ~ k h~ l pW {z" 1 ) , 
has the following expression^ 



(0.1) = M-} Z ~} S \ Y,T^ z n 

y ' V r n (t,s) 

2V ' V r n+1 (t,s) 

in terms of the 2n-uple integrals 

(0.2) r n (t,s):= [ [ dxdyA(x)A(y) f] e^T^' k - s ^+ cx 

J Jr 2 ™ IX 



kVk 



which form a vector r := (r n ) n >o of r-functions. This provides a concrete 
realization of the Sato representation of the two-Toda wave functions in terms 
of r-functions r n ; see [18], [2]. 

The pair of semi-infinite matrices L := (Li,L2)Jl defined by: 

(0.3) z^ x = z~ l % = Lj%, 



Ma] = («,« 2 /2,« 3 A...) 

2 Li is lower-triangular, except for a subdiagonal just above the diagonal with all entries = 1 
and L2 is upper- triangular, except for a nonzero subdiagonal just below the diagonal. 
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together with the pair (^1,^2)1 satisfy the two-Toda lattice equation; in par- 
ticular, 

r)T ■ r)T ■ 

(0.4) ^ = [(Lj) Uj Li] and ^- = [{I$)iM, i = 1, 2, 

for the usual Lie algebra splitting ( ) u , ( )i , explained in Section 2. The familiar 
standard Toda lattice (on tridiagonal matrices) is a special reduction of the 
two-Toda lattice, in the same way that the Korteweg-de Vries equation (KdV) 
is a reduction of the Kadomtsev-Petviashvili equation (KP). 

Conversely, starting from the the two-Toda lattice equations (0.4), one is 
led to wave functions and ^Z, an d a representation in terms of r-functions 
as in (0.1). As will be established in Section 3, the functions r n satisfy the 
standard KP-equation in t and s separately (see the beginning of §3), but they 
also satisfy another (new and useful) equation, which is third-order, relating 
t- and s-derivatives, namely: 

Theorem 0.1. Two-Toda r-functions T n (t,s) satisfy^ 

d 2 logr w d 2 log r n \ | d 2 logr n d 2 logr n | _ Q 
dt\ds2 ' dUdsi I I ds\dt2 ' dt\ds\ \ 

Vertex operators and "Christoffel-Darboux" kernels. Backlund-Darboux 
transformations refer to the general recipe of factorizing differential or dif- 
ference operators and flipping the factors, to form a new operator. When 
we let this situation flow in time, the new wave functions (eigenfunctions) 
can be expressed in terms of the old ones as Wronskians (continuous or dis- 
crete), and the new r- function is expressed in terms of the old ones, by 
means of vertex operators. So, the latter can be viewed as generators of 
Backlund-Darboux transformations for differential or difference operators at 
the level of r-functions. Typically, vertex operators X map r-functions into 
r-functions, and their squares vanish; although r satisfies a highly nonlin- 
ear equation, vertex operators have an additive property; that is r + Xr is a 
r-function as well! 

With the two-Toda lattice, we associate four different vertex operators 
Xjj(A,/i), for 1 < i,j < 2; they map infinite vectors of r-functions into r- 
vectors, as explained in Section 5. The vertex operators Xn and X22 are basic 
vertex operators for Toda, and KP, as well, whereas X12 and X21 are vertex 
operators, native to two-Toda. In particular, we construct 

X 12 Gu, A) = A- 1 x (A)X(-s, \)X(t, n)x(fi), 



3 in terms of the Wronskian {/, g}t = — fjjr ■ 
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with A the customary shift-operator (Av) n = t> n +i, and with 

X(t,X) :=eSi '* A V^ ~, , X (A) := diag(..., A" 1 , 1, A, . . .). 

Besides this work, the vertex operator X12 will also play a major role in our 
later work on symmetric and symplectic matrix integrals. Given a two-Toda 
lattice r-vector r = (...r_i, To, t\, . . .), we have that r + Xi2(y, z)r is another r- 
vector. But more is true. We show that the kernels K\ 2tn (y, z), defined by the 
ratios (Xi 2 T) n /T n , have eigenfunction expansions in terms of the eigenfunctions 
ty, reminiscent of the Christoffel-Darboux formula for orthogonal polynomials; 
to be precise, 

Theorem 0.2. We have (§5): 
(0.5) K l2yn (y,z):=-X 12 (y,z)T n = ]T ^(aT 1 )*^), 

Tfl -oo<j<n 

together with a Fredholm determinant-like formula, 

(0.6) det ((K ij>n {y a , z p )\< a «< fc = — ( f[ 3Qj(y/, z<>)t j . 

Tn \t=l ) n 

In the semi-infinite case, the sum in (0.5) is replaced by J2o<j< n - 

Vertex operators, Virasoro algebras and two- Toda symmetries. The vertex 
operators provide central extension realizations for Virasoro algebras; e.g., 
(0.7) 

^y k+1 X 12 (y,z) = [4 2) ,X 12 (y,*)] and JU fe+1 X 12 (y, z) = [if, X 12 (y, z)] 



with 

(0.8) Jf := 



'2n + k+ 1U (1) + n(n + 1U (0)> I 

nGZ 



= \ (jf + (2n + fc + 1)4" + n(n + 1) jf ) 



tl — »■ — S 



forming Virasoro algebras of central charge c = — 2. 

In (0.8), equals 5ko for £ = 0, Heisenberg generators for £ = 1 and 
Virasoro generators for £ = 2. It follows that the vertex operators Xjj of the 
two-Toda lattice form, upon expansion, the generators of a large algebra of 
symmetries, which come from the master symmetries for the pair of matrices 
L = (L 1 ,L 2 ). 

This is a special case (a = 1) of a more general statement concerning 
vector-vertex operators, depending on a parameter a, 

X a (t,u) = A~ l e£i ^e' a ^ ^^ x (n a ). 
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Theorem 0.3. The generators jj^(a), defined by 



dz 



z k+1 X a (t,z) 



^\a),X a (t,z) 



form a Virasoro algebra 
with central charge 

c= 1-6 



12 



a\ 1/2 fa\- 1 ^ 



Commutation of Virasoro and "vertex integral operators" . Consider now 
a more general weight p(y,z)dydz := pt, s (y, z)dydz := e Vt ' s ( y ' z ^dydz on R 2 , 
with po = e Vo , where 



oo oo 



(0.9) V t>a (y,z) := V (y,z) + £ t lV l - £ aj 

l l 

ij>l i>l i>l 

and a set -E C M 2 of the form, 

(0.10) E = E 1 xE 2 := U[ =1 [a 2i _i, a 2i ] x U| =1 [6 2i _i, 6 2i ] C K 2 , 

involving disjoint unions. The weight (0.9) and the boundary of the set (0.10) 
enable one to define two types of operators: 



(i) Virasoro- like operators, for k > — 1, 



d 







i=i 



with J^ 2) and j£ 2) as in (0.8). 



(ii) An operator defined by a vertex integral operator X12 over E, 



Ue := / / dxdy p(x,y)X 12 (x,y). 



(0.12) 

It is an important ingredient in this work that Vfc and (Ue)™ commute: 
[V fc , (U £ ) n ] = [V fe , (U £ ) n ] = for all n > 1 and k > -1. 
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Two-matrix integrals over product sets and Virasoro constraints. It follows 
that the 2ra-uple integral, like (0.2), but taken over the set E n = E'[ l x 
c R 2n , 

(0.13) T*(t,s):= [[ dfdyA(x)A(y)ffe^ ( ^^- s ^ )+ ^>i c -^, 
J jEn fc=i 

is related to r := (t? 2 ) , by 

V / n>0 

rf = ((U B ) n r) n . 
This implies, setting all Cjj = 0, but en = c: 

Theorem 0.4. r n and r,f satisfy the Virasoro-like partial differential 
equations, labeled by k = — 1, 0, 1, . . . : 

(0.14) (-j2a\ +l -^ + J { A^ + cp k ^)Pn{-d s )^ orti = 



i=l 



d 



db 



(2) ~(2) ~ ~ 

J k n and J), n were defined in (0.8), the Hirota symbols p{dt)q{—d s )fog in (6.1) 
and the p^s are the elementary Schur polynomials. 

Application to the spectrum of coupled random matrices. Consider a prod- 
uct ensemble (Mi, M 2 ) € : = 7Y n x7Y ra of n x n Hermitean matrices, equipped 
with a Gaussian probability measure, 

(0.15) C n dM 1 dM 2 e -^(^i 2 +M|-2cM 1 M 2)) 

where dM\dM 2 is Haar measure on the product 7Y^, with each dMj, 

n n 

(0.16) dMi = Al(x) Y[ dxdU and dM 2 = A 2 (y) ]J dy^C/ 

1 1 

decomposed into radial and angular parts. We define differential operators A k , 
B k of weight k, in terms of the coupling constant c, appearing in (0.15), and 
the boundary of the set 

(0.17) E = E 1 xE 2 := U[ =1 [a 2i _i, a 2i ] x Uf =1 [6 2 ;-i, M C M 2 . 

They form a closed Lie algebra, as spelled out in (11.4): 
(0.18) 
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The following theorem deals with the joint distribution (§11), 

(0.19) P n (E) := P(all(Mi-eigenvalues) G E 1 , ail(M 2 -eigenvalues) G E 2 ), 

and leads to a formula, which is the "mirror image" of Theorem 0.1. 

Theorem 0.5 (Gaussian probability). The statistics (0.19) satisfies 
the n- independent nonlinear third- order partial differential equations 
(E n :=l\ogP n {E)): 
(0.20) 

\B 2 AiF n , S^^+^-J^ - ^A 2 B r F n , AiB l F n + -^—^^ =0. 

Remark 1. Since the equation above for the joint statistics is independent 
of the size n, the same joint statistics for infinite coupled ensembles should 
presumably be given by the same partial differential equation. 

Remark 2. For E = E\ x E 2 := (— oo,a] x (—00,6], equation (0.20) takes 
on the following form: Upon introducing the new variables x := —a + cb, y : = 
— ac + b, the differential operators A\ and B\ take on the simple form A\ = 
d/dx, B\ = d/dy and (0.20) becomes 

d_ Ac 2 -l) 2 Sj + 2cx-(l+c 2 )^ d_ / V-l) 2 S + 2q/-(l + c 2 )x \ 

Remark 3. Equation (0.20) also has a "zero-curvature" formulation, namely: 
(0.21) [Ax-X n ,Bi-Y n ] = 0, 

with 

AiB±F n + and y "-^ 1 F„ + ^ T - 

The last section deals with coupled matrix ensembles, where the joint 
statistics is given by the "Laguerre distribution." Unlike Theorem 0.5, the 
Laguerre case for n x n matrices lead to (inductive) differential equations for 
the matrix integral (0.19); indeed, the equation contains a term, which is 
expressible in terms of the same expression for (n — 1) x (n — 1) matrices. 

Acknowledgment. We thank Taka Shiota for many useful discussions con- 
cerning Fay identities. We also thank Edward Frenkel for urging us to compute 
the central charge for the Virasoro algebra defined in Theorem 0.3. 



4 in terms of the Wronskian {/, g}x = Xf.g — f.Xg, with regard to a first order differential 
operator X. 
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1. Operators A and e with [A, e] = 1 and the 5-function 

Define the column vector \{ z ) = {z n )nez, and matrix operators A, A*, 
e, e* as follows: 

d 

^x(z) = zx(z), ex{z) = q^x(z), 



A*x(z) = z 1 x(z), e*x(z) = q^ZiX(z) 



Note that 



and 



(1.1) 



A* = A T = A" 1 , e* = -e T + A, 

A T x(z- 1 ) = z X (z~ r ), A X (z- V ) = z-'xiz- 1 ) 
e T X {z- 1 ) = z-'xiz- 1 ) - §- zX {z- 1 ) 
{ e^xiz- 1 ) = zxiz- 1 ) - e^rxiz- 1 ). 



(1.2) 



The operators A, A*,e, e* have the following matrix representation: 

A = (5ij-i)ij e z, e = diag(i) • A -1 = (i 5 ijj+1 )ij e z 
A* = (5ij +1 )i ij( zz, e* = -diag(z) • A = {-i8 i:j -i) i:je z- 
For future use, we also introduce the 5-function, 



(1.3) 

with the customary property 



1 t~ 

s(t)= y t n = + - 

w ^ 1-t 1-t- 1 



f(X^)S(-)=f(X,X)S(-: 

fl /i 



Note the function 
(1.4) 

has the usual ^-function property 

±-Jf(\ti)6(X,v)dti = f(\X) 
and is a function of A — fi only, since 

(1.5) < d 9 



For future use, we state: 
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Lemma 1.1 (T. Shiota). We have the following matrix representation^ 
(1.6) e^- A ) £ <5(A,A) = ~x(m)®X*(A) and 



A 



Proof. Note that, sincei 

A™ = (5ij- n )ij & Z, £ n = ((i)n<5i,j+n)i,jeZ , 

A* n = (5ij +n )i : j & z, £* n = ((—i)n^i,j-n)i,jeZ , 

we have 

5(A, A) = J] A^A-"- 1 = ^ A™(4 j+n+1 ) UeZ = (A^" 1 )^ 

<5(A, A*) = J] A"A" +1 = A n (<5v-n-i)ijez = (A^ -1 )^ , 

and 
Hence 

e^- A)£ 5(A,A) = | £ ( _^ k ) ((j, - A)'- fc A*- J '- 1 

i-fc>0 



5 Given two column vectors a and 6, the matrix a <g) b is defined componentwise as follows 

(a (g> 6)y = Oj6j. 

6 For any k,n £ Z, n > 0, we use the standard notation (k) n := — l)...(fc — n + 1) and 
fc I _ (*)« 
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2. The two-Toda lattice 



Consider the splitting of the algebra T> of pairs (Pi, P 2 ) of infinite (Z x Z) 
matrices such that (P\)ij = for j — i » and (P2)ij = for i — j » 0, used 
in [6]; to wit: 

V = V + + V_, 
V+ = {(P, P) j Pij = if \i - j\ » 0} = {(Pi, P 2 ) € V J Pi = P 2 }, 

25- = {(A, P 2 ) I (Pi)ij = if j > i, (P 2 ) ij = if i > j}, 
with (Pi,P 2 ) = (Pi,P 2 )+ + (Pi,P 2 )- given by 

(2-1) (A,P 2 )+ = (Plu + P2t,Pl» + P2t), 

(Pi,P 2 )_ = (Pu-P2e,P2u-Piu); 

P u and P^ denote the upper (including diagonal) and strictly lower triangular 
parts of the matrix P, respectively. 

Throughout this paper, we will use the following operators (a mul- 

tiplication operator) and e Vi ^ (a shift), where 



(2.2) 



1 

00 d 

*(*) = E — ar 



and 
and 



6W = E^ Z fc ' 



00 3 



so that 



e °^ +fc »/(t, a ) = /(t + a[z- 1 ], a + 6[z]) 
with [a] = (a, a 2 /2, a 3 /3, . . .). 

The two-dimensional Toda lattice equations 



(2.3) 



dL 

dt n 



(^,o)_ 



and 



dL 

ds n 



(o,^) + , 



n = l,2,... 



are deformations of a pair of infinite matrices 
(2.4) L = (L 1 ,L 2 ) = ( J2 ^A*, E "f'A')^. 



-oo<i<l 



-l<i<oo 



with A the shift operator of Section 1 and where aj 1 ^ and are diagonal 
matrices depending on t = (ti, i 2 , . . .) and s = (si, s 2 , . . .), such that 

aS 1} = / and (a^ ^ for all n. 

V / nn 

In analogy with Sato's theory, in [18] it is shown that a solution L of (2.3) 
has the representation 

Li = WiAV^f 1 = SiASf 1 , L 2 = W^A -1 ^ -1 = 5 2 A -1 5^ 1 
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in terms of two pairs of wave operators 

f S t = E,< Ci(t, s)A\ S 2 = E,> <(t, s)A i 

[ Cj,c£ : diagonal matrices, cq = I, (c' )u ^ 0, for all i 

and 

(2.5) W l = S i (t,s)e^ A \ 

One also introduces pairs of wave and adjoint wave vectors ^ = (^i, ^2)1 and 

(2.6) ^i(t,s;z) = W iX (z) = e^S iX (z), 
*Ut,s;z) = (W?)- I x *(z) = e-^ Z \SJ)- I x*(z), 

which evolve in t and s according to the following differential equations S 



(2.7) 



^* = (L?,0) + * = ((L?) U) (L5f) w )* 

^* = -((L?,0) + ) T vr 
J-** = -((0,L2) + ) T **. 



Besides L = (Li,L 2 ), we define the operators L* = (L*,^), M = (Mi,M 2 ) 
and M* = (M{,M%) as follows 

(2.8) 

L := (W1AW1 1 , W 2 A*W 2 ~ 1 ), L* := ((Wi)~ 1 A*W?, (W^^AWj), 
M := (WxeWi \ W 2 e*W^ 1 ), M* := ((W^)~ 1 £*W^, (W 2 T )- 1 eW 2 T ), 

which satisfy, in view of (2.6) and (1.1): 

L* = (z,z- 1 )*, Mtf = (j-^^y, [L,M] = (1,1), 

L*V* = (z,z- 1 )**, M*tf = g^W, [L*,M*] = (1,1). 



7 Here the action is viewed componentwise, e.g., (A, B)ty = (A^i , B^2) or (2,2 1 )>I' 
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The operators L, M, L T , M T and W := (W±, W 2 ) evolve according to 



(2.9) 



JL L 
dt n I M 



a 
dW 



L 

M 



(0,^5)+, 



L 

M 

L 

M 



(L?,0)+W, and 



dW 



dt n K 1: ^ ' ds n 

Ueno and Takasaki [18] show that the two-Toda deformations of and 
hence L, can ultimately all be expressed in terms of one sequence of r-functions 

r(n,t,s) = T n (t 1 ,t 2 ,...;s 1 ,s 2 ,...) = det[(Si 1 S 2 {t,s)) i j]- 00 < i j< n - 1 , neZ: 

to wit: 



(2.10) 



*i(t,s;z) = 

V 2 (t,8;z) = 

n(t,s;z) = 

V 2 (t,s;z) = 



T n (t,s) 



e '«T„+i(t, 



e r * Tn+1 (t,s) c _ T ~ t . 

T n+1 (t, s) 



e m T n (t, s) _y-°c 

e ^— 'i 

r n+ i(t, s) 



Finally the pair of matrices = (Wi,W"2) satisfies the bilinear relation 
(in the ± splitting of (2.1)) 

(W(t,s)W(t',s')- 1 )- = 

or equivalently, 

(2.11) Wi(t, s)Wi(t', s')" 1 = W 2 (t, s)W 2 (t', s')-\ 

from which one proves Proposition 2.1; for details see [6]. Equation (2.13) 
below is established in Adler-van Moerbeke [4]. 

Proposition 2.1 (bi-infinite and semi- infinite). The wave and adjoint 
wave functions satisfy, for all m,n G Z (bi-infinite) and m,n > (semi- 
infinite) and t, s, t', s' G C°°: 

(2.12) 



./z=oo 2/KIZ J z =0 



I J.js dz 



2iriz 
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Moreover L\ has the following representation in terms of r-functions 

(2.13) L\ = V diag ( P^/dt)r n+k . e+1 or n \ ^ 

E=0 \ T n+k -e + iT n J n& 

In the semi-infinite case, we reinterpret A -1 as A T , where A is the semi- 
infinite shift operator. Then one shows (2.11) and (2.12) are also valid. Also 
the semi-infinite case is obtained from the infinite case by setting r_j = for 
i > 1 and To = 1. Then the semi-infinite wave vectors 

(*in(i, a; *)e" and (% n {t, s; z)e^ 

V / n>0 \ / n>0 

are vectors of polynomials of degree n = 0, 1, 2, . . . in z and z~ l respectively, 
as follows from (2.12); see [2]. In the semi-infinite case, we must define L\ and 
Lj n for integers n < 0; namely for neZ: 

(2.14) 4 n) (e - £ *****!(*)) = TT+C^e-E^tfiC*)), 

where 7r + refers to the projection ft+(J2iez a i z% ) = Si>o a « zi - Observe that 
= L™, L% = Lj™ for n G Z, n > 0; indeed, multiplying the vector 
^i(z)e _s * i2:l of polynomials with z n , n > maintains, the polynomial charac- 
ter, and thus for n > 0, 

LS n) (^!(z)) = 7r+(*"*i(*)) = z^i(z) = L?tt(s). 



3. Bilinear Fay identities and a new identity 
for two-Toda r-functions 

Two-Toda r-functions r(i, s) satisfy the KP-hierarchy in t and s sepa- 
rately, of which the first equation reads: 

o \\ (fd\ 2 , V / s \ 2 d 2 



wj iogr+ %u; iogr j +3 uj iogr - 4 ^ iogr=o - 

But they also satisfy the following identity: 

Theorem 3.1. Two-Toda r-functions satisfy^ 

(<9 2 logr n d 2 logr n l [<9 2 logr n <9 2 logr n l 



(3.1 



dt\ds2 ' dtidsi [ ] ds\dt2 ' dt\ds\ \ 

) t\ \ - ) si 



8 in terms of the Wronskian {/, g}t = — / §f ■ 
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The proof of this theorem hinges on the bilinear identity, due to Ueno- 
Takasaki [18] and a number of lemmas: 

Proposition 3.2. Two-Toda T-functions satisfy the following bilinear 
identities: 

(3.2) / r n (t - [z-\s)r m+1 {t' + [z~\ s')e^ (U-^ z n-m-l dz 

J 2 = 00 

r il+1 (t,s- [z})r m (t',s' + [z^^-^z^-'dz, 



lz=0 

or, expressed in terms of the Hirota symbol^ 

(3.3) £p m _ n+ ,(-2a)p,(^)e E " K ^ +bft ^ ) r m+1 o r n 
3=0 

= Y,P-m + n +j (-2b) Pj (d s )e^ (afe ^ +bfc ^ > Tm o Tn+1 , 
3=0 

both, for the bi-infinite (n,m G Z) and the semi-infinite case [n,m G Z, 
n,m> 0). 

Proof. (3.2) follows at once from Proposition 2.1 and the r-function rep- 
resentations (2.10), whereas (3.3) follows from the shifts t i— > t — a, t' \— * 
t' + a, s h s - i, s' i— > s' + b, combined with the definition of the Hirota 
symbol. □ 

Proposition 3.3. 

,5 \ d 2 logTn + l 

/ 34 n / L fc\ _ Pk-l(Ot)T n +2 °T n = dsidtk 

V Vn.n+l r n+2 r n 82 '°S!>+ 1 

/ 5 \ 9 2 logr n+ i 

(3.5) (/^/r 1 ) = Pfc - l( " as)r " +2 ° Tn = W^^. 

V )n,n+X T n+2 T n 9 i° g I" +1 

Proof. Set m = n + 1, all 6^ and = 0, except for one aj+i, in the Hirota 
bilinear relation (3.3). The first nonzero term in the sum on the left-hand side 
of that relation, which is also the only one containing flj+i linearly, reads 
(3.6) 

p j+1 (-2a)pj(d t )e 3+1 9t ^+ 1 r n+2 o r„ + ... = -2a j+ iPj(dt)T n+2 or n + 0(a 2 j+l ), 

whereas the right-hand side equals 

a 9 d d 
(3.7) p (p)pi(d s )e J+lat ^+ 1 T n+1 oT n+1 = — (l + a j+ i— \- ...)r n+1 o Tn+1 . 



d Sl y ' J ^dt 



'For the customary Hirota symbol p(dt)f o g := p(J^)f(t + y)g(t — y) 



y=0 
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Comparing the coefficients of Oj+i in (3.4) and (3.5) yields 

d 2 

-2 Pj(d t )T n+2 °T n = Q g ^g t - Tra+1 ° T n +r, 

in particular, we have 



Pk-l(dt)T n +2 °T n _ 2 

and so, for k = 1, 



3.8) 2 = ~ a ^gT n+1 , 

t 2 +1 dsxdtk 



(3.9) TnTn+2 = _ " logT„ + i- 

T n+l ds l dt l 

Dividing (3.8) and (3.9) leads to the second equality in (3.4). But, according 
to (2.13), the (n, n + l)-entry of L\ is given by (3.4). The similar result for L\ 
is given by the involution 

t < — > —s and L\ < — > hLjh^ 1 . □ 

Proof of Theorem 3.1. Set k = 2 in the identities of Proposition 3.3; then 
subtracting of identity (3.5) from of (3.4) leads to Theorem 3.1. □ 

The (n,n + l)-entries of L\ and /iLj 2 /i _1 have the following equivalent 
expressions, which will be useful in the theory of Toeplitz matrices, as applied 
to the distribution of the length of the longest increasing sequences of a random 
permutation. The second identity, appearing in L 2 below is an expression 
purely in terms of one component T n+ ±, at the expense of introducing a djdti- 
derivative; the third identity involves d/dt\ and d/ds\ only, but at the expense 
of involving nearest neighbors r n and T n+ ±. 

Lemma 3.4. Two-Toda r-functions satisfy: 

9 i„„ T n+2 



d 2 



(3.10) (L?) = ^log 

V / n,n+l Oil T 5 



n 



(hLl 2 h^) 



n,n+l 




log r n+ i 




5 Sl ^ l0gr " +1 
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Proof. From Proposition 3.3 (k = 2), it follows that 

d 2 logr n+ i d 2 logr ra+ i |(%2°T n ) 

dsidt 2 dsidh T n+2 T n 

d 2 logr„ + i d r n+2 
ds\dti dti r n 

a 2 logr n+ i / d T n+2 d r n+ i 



— log 1 log ■ 

ds\dt\ \dti T n+ \ dti r, 

= d 2 logr ra+ i / d / r ra+ i d 2 ^ \ 
ds\dt\ \dt\ I r n ds\dt\ n + l J 

+ 4r lo g , using (3.9) 

d 2 \ogr n+l ( d r n+1 d ( d 2 \\ 

9 a r n+1 d 2 d ( d 2 \ 

ati r„ osioti oti yosioti J 

The second to the last equation establishes the first equation (3.10). The 
second equation (3.10) is simply the dual of the first one by ti <-> — s,-. □ 

The remaining statements in this section hold for both the bi-infinite case 
and the semi-infinite case; we thank T. Shiota for showing us how shifting 
the arguments in various directions, and repeatedly, leads to many different 
identities. 

Proposition 3.5 (the Fay identity). // 

(3.11) a-a' = f:^ 1 ] - E^fc 1 ] and P-P' = EM " !>*]> 
l 

with p,q,p',q' > 0, we have 



11 11 



(3.12) f> n (a - [z;\{3)T m+1 {a> + \zj\0)z\^ "f= l( ^_i _ 

^=i llfc=il^ ^fc J 



^(r„( Q -H,/3)r m+l(a ' + H,/3')|| r |^ 



x=0 
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-r>-illl = i(ve-u k ) 



U P k=i {n - v k ) 



r'\ dx r 



Eli (x - v k ) , 



x=0 



where r:=n — m + q— p and r' := — n + m + q' — p', with the understanding 
that (d/dx) r = for r < 0. 



Proof. The relations (3.11) imply 



(3.13) 



n?(i-£) 



and e*™-®*- 



n| (i - j) 
rrf (i - ^) 



If / denotes a holomorphic function in a large enough disc around z = oo, 
as in 

f(z) := r n (a - [z' 1 ], P)r m+1 (a' + [z' 1 }^'), 
then we have, using (3.13), 

I f{z)e T '^^-< )zl z n - m ' 1 dz 

Z7TI J z=oo 

i r ,n?fi 



2vri 



n?(i-i) 



/ £n— m— 1 



-n+rn- 



-q+p upon setting 



= 1 



rifc=i (^ 1 z k 1 ) 



r! \dx 



x=0 



which uses the fact that the integrand has poles at x = z^ 1 (1 < k < p) and 
at x = 0, if r = n — m + q— p > 0. In = the sign change in dx due to x = 1/z 
and the change of orientation of the contour integration cancel each other out. 

If /' denotes a holomorphic function in a large enough disc around z = 0, 
as in 

f'{z) := r n+1 (a,(3 - [z])r m {a' , + \z\\ 



then we have, again using (3.13), 
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— / f{ Z )e^-^- i z n-m-l dz 

2m J z= o 



r\i (i _ y±\ 
*=° U{ (1-^) 



1 

2^i 

2mf z=0 J >Ytf(z-v k ) 

Ew/ , >, n|=lC^ ~ "fc) n-m+p'-q'-l 
<=i rifc=i - «*:) 



i 



rii (a;- 1 "*) 



Eli (a - «fc) 



z=0 



which used the fact that the integrand has poles at z = v k (1 < k < p') and at 
z = 0, if r' = —n + m — p' + q' > 0. □ 

Corollary 3.6. 

T n (t - [Z^ 1 ], S + [v] ~ [u])T n (t, S) - T n (t, S + [v]~ [u])T n (t - [z~\ s) 



V — U 



T n +l(t, S - [u])T n -l(t - [Z ],S + [v]). 



Proof. Setting m = n — 1, a = t, f3 = s + [v] — [u], a' = t — [z -1 ], (3' = s, 
we have a — a' = [z^ 1 ] and — f3' = [v] — [u], and thus p = p' = q' = 1, 
q = 0, with = —p' + q' < n — m = p — q = 1; that is, r = 0, r' = —1. Then 
Proposition 3.5 leads to the proof of Corollary 3.6. □ 

Corollary 3.7. 

z -i 

T n (t,S+ [vi])T n+1 (t + [Z X X ] -[Z 2 \s- [f 2 ])— T i — 

z l ~ z 2 

+ r n (t+ [z ± x ] - [z 2 1 ],s + [vi])T n+1 (t, s - [v 2 })^y^ — ZT 

z 2 — z l 

Vl 



T n +l(t + [Z X s)T n {t -[Z 2 1 ],S+ [Vl] ~ [V 2 ])- 



vi - v 2 

v 2 



+ T n+1 (t + fa 1 ], S + [ui] - [v 2 ])T n (t ~ [Z 2 1 ],S) 

v 2 — V\ 



Proof. Setting m = n, a = t + [z± ], a' = t — [z 2 ], (3 = s + [v\], (3' = 
s — [v 2 ], we have a — a' = [zf 1 ] + [z 2 l \ and (3 — (3' = [v{\ + [v 2 ], and thus 
p = 2, q = 0, p' = 2, g' = 0, with —2 = —p' + q'<0 = n — m<p — q = 2, and 
so, r = r' = —2. Similarly Proposition 3.5 ends the proof of Corollary 3.7. □ 
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4. Higher Fay identities for the two-Toda lattice 



Lemma 4.1. 



fc+i 



E T N -i(t -[z f 1 ],s + [yJ])T N - k 



( fe+i k+i \ 
f-Et^WEfor 1 ] 



v 



J 



(k+l k+l \ fc 

* - E fe -1 ]. s + E br 1 ] n (y™ - y'i 
i i / <=i 



Proof. In Proposition 3.5, we set 
a = t, 



fc+i 

j'=i 
fe+i 



i=l 



Now obviously, 



fe+i fc+i 
a - a' = Et^ 1 ] /?-/?' = [y,; 1 ] - E fof 
i=i 



i=l 



p = A; + 1, q = 0,p' = 1, q' = k, n = N — 1, m = N — k — 1, r = r' 
With these data, we have 



k+l 

E 

e=i 



E^v-i(«- 1 ],/3)r JV _ fe (a / + [z e \(3') . _ 1 — 



fc+i 



= r^a, - [ J /- 1 ])r JV _ fc _i(o', + [y- 1 ]) JJ (l/m " i/f *) 
establishing Lemma 4.1. 
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Theorem 4.2. The two-Toda tau-functions r(i, s) satisfy the following 
two higher Fay identities: 

'TNit-iz-^ + iy- 1 }^) 1_ 

Uj — ZiJ i<i,j<k 



(4.1) drtf 7 *"-*' 4 ^ 1 ) 

\ T N (t,S) Vj-ZiJ 

,M-I> A(y)A(z) ^(* + Sfb/r 1 ]-E^7 1 ],- 



= (-1)- 



Ilfc,/(f* -2 *) rjv(t,s) 

/ r/v_i ft — b. -1 !. s 

det 



Vjv_i(t- [zr\ s + [ y] Y 



Tjv(M) A<ij<fc 
Aiy-^Aiz- 1 ) ^ 



T N (t,s) 

Proof. The inductive method for proving the first identity is due to [6]. 
As to the second relation, we also proceed by induction on the index k. Since 
the identity is obviously true for k = 1, we assume it to be valid for k > 1 and 
we prove its validity for k + 1. Indeed, by expanding the determinant according 
to the first column, we find 



det 



TAT-ift-fcrV + ty, 71 ]) 

T N (t, s) 



l<i,j<k+l 

~ ^ } r N (t,s) de M r N (t,s) 



= i 



fe+i 



E, ,,e-i TN-i(t - [z t 1 },s + [y l x ]) TT / -l -k TT , -i -is 

^=1 i<i<j<fc+i i<i<j<fc+i 



T N {t, S 

Aiy-^Aiz- 1 ) 



, by induction 



2J '•' 



= A^-^A^" 1 ) ^ — 

tjv(i, s) 

by Lemma 4.1, ending the proof of Theorem 4.2. □ 

Such Fay identities were also obtained in the context of the multicompo- 
nent KP hierarchy by J. van de Leur [19]. 
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5. Eigenfunction expansions and Vertex operators 

In terms of the vertex operator, 

x(t,\) := e sr^A- e -Er A " 



acting on functions /(ti,t 2! ---) of t € C°° and using the diagonal matrix 
x(A), define the following four operators acting on column vectors g = 
{.9n(h,t2, ■ ■ -))nez : 

(5.1) Xi(t, /i) := X(t, fi)x(fi), Xj(t, A) := - X *(A)X(-i, A), 

X 2 (,s, M ) := -X(s, m)x*(m)A, X$( S , A) := A'V W"*. A), 



and the compositions 
(5.2) 

/ Xn(/x, A) X 2 i(/i,A) 
1 X 12 (//,A) X 22 (/x,A) 



The main theorem of this section is the following 

Theorem 5.1. TTie following holds: (i) in £/ie bi-infinite case, 



(XJ(t,A) Xl( S ,A))®(X 1 (t,/x) X 2 (s,/x)) 

/ XJ(*,A)Xi(t,/i) Xf(i,A)X 2 ( S , M ) \ 
^ X?( a> A)Xi(t,M) X|( S ,A)X 2 ( S , M ) J ' 

Jio 



(5.3) 



(E i <n^ i (A- 1 )^ lj (M))n e Z (E i >n*2 i (A- 1 )*2 i (M- 1 ))n6Z 
= 1 / Xn(/X,A) X 21 (/X,A) ^ 

r I Xi 2 (/i,A) X 22 (/x, A) 



and (ii) in the semi-infinite case 
(5.4) 



(Eo<i<n *y(A)*y (M))n>o (£,•>„ ^(A)^/^ 1 ))^ 

(Eo<,<n*2i(A- 1 )*li(M))n>0 (Ei^t^A" 1 )*^^- 1 ))^ 



1/fXii X 21 \ / (i-M)-i e S^(M-^ 
r 1 1 X12 X 22 J + I 



1 For column vectors t>i, 112, uii and 102, we define 
(vi v 2 ) ® (uii iu 2 ) = 



iii 5< li'i 5s 11)2 

l>2 ® «Jl W2 ® ?"2 
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Corollary 5.2 (eigenfunction expansion). In the semi-infinite case, the 
functions below admit the following eigenfunction expansions: 

oo 

(5.5) (1 - ^)-i e Sr*™(M"-A") = £*J.(A)vl/y(M) H<|A| 

A 3=0 

(i - V^r^-A-) = E^^ 1 )^^- 1 ) 

^ j=o 

oo 

i=o 

The proof of Theorem 5.1 and Corollary 5.2 relies on Lemmas 5.3 and 5.4. 
Remark 1. With the identities 

oo 

(5.6) g-Er^ = l- a and ^ a* = (1 - a)" 1 , 

o 

the composition of X(t, fi) and X(—t, A) relates to the customary vertex oper- 
ator X(t, A, (jl), as follows: 

(5.7) X(-t,X)X(t,f,)f = X(-t,\)(e S ^ i f(t-[ f ,- 1 ])) 

e a l (^-v) /(t+[A -i ] _ ^-1]-, 
N X(t,»,\)f(t), 

A — /i 

where X(t,ji, A) admits the following expansion in terms of VK-generators: 

(5.8) *(t jM ,A) := exp^t i (^-A l )jexp^(A- l - / .- l ) 7 — 

OO f \ \ oo 

= E A-^ (fc) , with Wf= fo. 

fc=0 ' £=-oo 

Note that In and X22 are closely related to the two-Toda vertex operators 
defined in [6]. Acting on infinite vectors of r-functions,0 they give 

(5.9) X U ( M ,A) =XJXi = - X *(\)x(v)X(-t,\)X(t,ti) 



e 

A 

A — fj, 
A 
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^ — A V V A/ / n( zz 



^ A Vfc=0 f=-oo /nGZ 



-X(i,A,/i), 



/x — A 
(5.10) 

X 22 ( M ,A)=X*X 2 = -A- 1 x (A)x*(/i)X(- S ,A)X( S , M )A 



'' X(.,/x,A) 



/x — A \ \ /x 



/ OO / \ \ oc N 



-X(s, A, fi) 



k=0 ' ' l=-oo / nGZ 

/J 



wit! 



12 



/x — A 



<i = E ( • V kWt j) and W«=WW 



j=0 



J 



Remark 2. One easily computes from (5.8) and (5.9): 
WP = S nt0 , W^ l) = 4 X) and = J® - (n + n G Z 

<9/<% n if n > 



. r (i). 

if n = i+i= " 



jP-.= { (-n)t_ n ifn<0 , 4 2 ):= £ rJ^J, 



and 
(5.11) 

W^ 1 ] =wl l) +mW^ ] W%1 =wl 2) + 2mW^ l) +m{m-l)wf )) 

= + mSio , = 4 2) + (2m - * - 1) jf } + m(m - l)fc . 

Before establishing Theorem 5.1 we first prove the following lemmas: 



2 Note that in the notation of [6] 



£(i,A,A*) = - — -X u (/i,A) and X(s,A,/i) = - — -X 22 0, A). 
A fi 
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Lemma 5.3. The following hold: 

(i) *!,„(A)* 1 ,„( At ) = - Hz!^^ + (£)" *(- t ,,)x (t , M W 

(ii) *2,n(A _1 )*2,n(/i _1 ) = (£)" ^£^1^ _ *(^)*(-..A)t„ +1 

(iii) *2,n(^ _1 )*l,n(A t ) = (M)" ^M^iilln _ (/xA) n-l j^A^fOTn^ ^ 

(iv) *i,„_ 1 (A)* 2 ,„_l( A .- 1 ) = (mA)- +1 H^Hif^ _ (/tA) -» Xi-t^XMr^ 



Proof, (i) By the explicit expression for X(t, A), the right-hand side of the 
first relation equals 



(ay 1 -srw-A*) f / 



A T„+l(t, S) T„(t, S) 



(^\ n c T,ru^-^) r n+ i(t+[\ \s)T n {t-\n 1 ],s) 

\XJ T n+ l(t,s)T n (t,s) 

= *J,„(A)tti,„(/i) 

by Corollary 3.7 with t>2 — ► 0, — ► 0, zi = A and Z2 = /x, and finally by 
(2.10). 

(ii) Similarly 

A" gSr^'-A') 1 f A rn+ift.s+tA- 1 ]-^- 1 ]) | ^(t.s + tA- 1 ]-^- 1 ]) 
A*/ ' 1-^ \ /" r n+ i(t,s) r„(i,s) 

A \ " gEfsiC^-A*) T n (t,s + [A~ 1 ])r» + i(^,s — [p" 1 ]) 

A* / T„ + l(t, s)T„(t, S) 

= *;,n(A" 1 )*2,n(^ 1 ) 

by Corollary 3.7, with z\ — > oo, Z2 — > oo, t>i = A -1 and t>2 = and by 

(2.10). 

(iii) Also, the right-hand side of the third relation equals 

M » e Er(^"-M') ( Mt-lfi-'ls+lX- 1 }) _ 1 r w -i(t- [M'^g+fA- 1 ]) 

\ Tn+l(t,s) ^iA T„(t, s) 

= ( M A)"c E? ° (t ' M '~ s ' A ' ) ' r " ( ' f ' S + [ A ~ 1 ]) Tn ( t ~ [m' 1 ]^) 

r„+i(t, s)r„(i, s) 

= *2,n(A _1 )*l,n(A«) 

by Corollary 3.6 with u = 0, v = A -1 and z = /j,, and by (2.10). 
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(iv) Finally, 

( A(A )-»+i c £rW-t«A i ) ( 1 r n+1 (t+[A- 1 ], S -[ M - 1 ]) + r^t+lX-^s-lfi- 1 ]) 

\ fj,\ T„(t,s) T n -l(t,s) 

(uA)~" +1 c E? ° (s ' M '~ t ' V) r "^ + [ A ~ 1 ]' s ) r "( t ' s Z t^' 1 !) 

T n (t,s)r„_i(t,s) 

= *t, n _ 1 (A)* 2 ,„-i(M" 1 ) 
by Corollary 3.6 with t i-> i+[A _1 ], z = A, v = and u = /i" 1 and by (2.10). □ 

In the next lemma we show that the Christoffel-Darboux type kernels, 
formed by means of the two-Toda wave function (2.10) can be expressed in 
terms of vertex operator acting on the r-functions; set X := X(/x, A): 

Lemma 5.4. When 1 < |/iA| < \fi\ 2 , the following holds (bi-inftnite case) 



(i) 




-(£)' 


(ii) 


E^n^A- 1 )^^- 1 ) = 




(iii) 






(iv) 




(ma)- 



i = (r-^ia^Jn, 



The proof is based on summing up the expressions in Lemma 5.3 and 
noting that, given the inequalities above, 



( — 1 and (/uA) n — > when n — > — oo 



-A- 
and 

/ X \ n 

and (//A) n — > when n — > +oo. 



A 
// 

Proof of Theorem 5.1 and Corollary 5.2. In the bi-infinite case, the state- 
ment of Lemma 5.4 leads at once to (5.3), whereas in the semi-infinite case, 
summing up the expression (i) of Lemma 5.3 yields a boundary term, by the 
fact that ro = 1. For (iii) the boundary term vanishes, since r_i = 0. To 
prove the first expansion of the corollary, let n | oo in the (l,l)-entry of (5.4), 
assuming |///A| < 1. Setting n = in the (2,2) and (l,2)-entries of (5.4) yields 
the second and third relations of (5.5), after first stripping the exponential part 
from the equation. □ 
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6. A remarkable trace formula 

Given two differential polynomials pipt) and q(d s ), define the customary 
Hirota operation: 

d d 

(6.1) p(dt)q(d s )fog:=p( — )q( — )f(t + y,s + z)g(t-y,s-z) 

dy oz y=z=0 

Theorem 6.1. In the semi-infinite case, we have the following trace 
formula involving elementary Schur polynomials, for n,m > —1: 

(6.2) L 2 + % = -4-^Pn+N(d t )p m+N (-d s )n O rjv-l. 
Remark 1. Here is an alternative way of writing (6.2): 

(6.3) 

]T (L? +1 1%+% = -4r^ E MPi>(dt)pA-ds)T N -i(t,s) 
0<i<JV--i TN ^ S > 

— — j+j =m+N 

where the the moments 

IHj ■= J J x i y^e v( - x ^dxdy. 
We shall need the following matrix of matrix operators:^! 

(6.4) 

( Nu N12 \ _ ( W 1 e^- x ^ £ 5(A/X)W{' 1 Wie^ 1 ~^ e * <5(A* / A) W^ 1 \ 
V N 2 i N 22 ) ~ \ -W 2 e { ^ 1 -^ £ 5{K/\)W{ 1 -W 2 e^- x ^*5{A*/X)W 2 - 1 ) 

where 5(A/X) = XS(X, A), with the delta function S(z) = Y^w.=-oo zTl i as defined 
in (1.3). 

Proposition 6.2. The following holds: 

(6.5) 

( Nu N l2 \ ( *i(/i)®*I(A) *i(a»)®*5(> -1 ) \ 

V ^2i ^22 y ' -*2(m -1 )®**(a) -^(m -1 )®*^* -1 ) ; 

Proof. Using (1.4) and Lemma 1.1, we compute 
Nn = XWie^~ x)e 8{\A)W^ 

= {W xX {p)) ® ((Wi T ) -1 X*(A)) = *i(M) ® *I(A) 



13 Note that in the notation of [6], N\ = if T ^A r n and jV 2 = - e -An 2 2 
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and 

N 22 = -\W 2 e ( ^- A)e *5(A,A*) W 2 l 

= -w 2X *(m)®x(a) w 2 x 

= -(WaxCA*- 1 )) ® ((W 2 T )- VtA" 1 )) = -^(a*" 1 ) ® ^(A- 1 ). 
The remaining relations are established in a similar way. □ 

Remark. The operators N(t,fi,X) have been considered in [6]; they are 
generating functions of symmetries on the ^-manifold in the following sense, 
by (2.8): 

" fe=l ' ^=— oo 

Proposition 6.3. In i/ie bi-infinite case, the following holds: 

(6.6) Yl z ~ nL i = 

5>- n L£ = * 2 (^- 1 )®*5(^- 1 ). 

Proof. Setting // = A = z in formula (6.4) for Nu, we find 

iV«(t, z, z) = {-ly-HiU/z) = {-If' 1 ]T z~ n Lf, 

which combined with Proposition 6.2 yields Proposition 6.3. □ 

Proposition 6.4. In the semi-infinite case, the following holds (see 
notation (2.14)): 

(6.7) 5>-4 n) = 

Y^ n L 2 (n) = *5(«" 1 )®*2(«- 1 ). 
Proof. Acting with J2 n 

& z^L^ on the wave function ^i(z), recalling 
that 7r + is the projection 7r+(J2iez) = J2i>o o-i^i an d using the usual property 
of the (5-function in the third equality, we obtain: 

(E*r4 B) W*)e- E " ti * < = E^ + (^i(^~ Er ^) 
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Zl 



^Oe-Er^ E(- 



1 

21 



Eoc , i P^Jl ti( zl 2 l) 



Z 
Z\ 



j=0 

by (5.5) 



which is valid for all z, z\ £ C; this establishes the first relation of Proposi- 
tion 6.4. 

Similarly, one shows for all z, Z2 6C 



z/z 2 



j=0 

by (5.5) 

= m(z^) ® v& 2 (^ 1 ))^(z- 1 ) e sr^ ) 

leading to the second relation. □ 

Proof of Theorem 6.1. Multiplying the first relation of Proposition 6.3 
with the second transposed, and using Proposition 6.4 in the second equality, 
we find 

£ A~>-™4 n) (Lj (m) ) T = ^ A-"4 n) £ ^ m (4 (m) ) T 

= (^i(A)^^t(A))(^ 2 (^ 1 )^^(^ 1 )) 
= *i(A) ® *2(A t_1 )(*i(A), *2(At -1 )), 



using regular matrix multiplication, in the last equality. 
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Upon taking the trace of the matrix above up to TV — 1, and using both 
Theorem 5.1 and Corollary 5.2, we find 

(6.8) r N (t,s) £ A-V™ E (4 n) (4 M ) T k 

n,m£Z 0<i<N-l 

= T N (t, S) E *li W*2i(A* _1 ) • <*l(A), *2GO) 
0<i<AT-l 

= (A^-^^-w")^^ - [A" 1 ], s + f^ 1 ]) 
. gEJW-fnA")^ + [A-i], s _ [^-i]) 

= (A/xJ^-ViCt + [A" 1 ], s - [/x-^riv-iCi - [A- 1 ], s + f^ 1 ]). 
Then, using the Taylor expansion, in A" 1 and /x -1 , we find: 



n=0 

oo / oo \ 

= E p n (±5 t ) E Pm(=F sV~ m A"" 

n=0 \m=0 / 

= f] ( Pll (±9 t )p m (T9 s )/(M))A->- m . 



m,n=0 

Therefore, on the one hand, 

oo 

(6.9) Tl (t + [A- 1 ], s -[ / x- 1 ])= ]T P«(ftW(-4)ri(M)A"' l /i- m ; 

m,n=0 

on the other hand, using the explicit matrix representation of t\ (see [2]), we 
see that 

(6.10) 

Tl{t+[X -\ s -[^]) = J J e ^ T{U+ ^ )xl - {s ^ )yl) e v ^dxdy 



= II (T^fKU-J) * iv 

J2 A~>~ m f ( x n y m e Vt ^ x ' y Uxdy 



m,n>0 
m,n>0 

which, upon comparison with (6.9), leads to 

Hn,rn{t, s) = p n (B t )p m {-& s )ti (t , s). 
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Therefore (6.8) reads 

TN (t,s) y, ^>' rn E (4 B) (^ (m) ) T )« 

n,meZ 0<i<iV-l 

= (Xfif- 1 YPi(dt) Pj (-d s )T 1 (t,s)\- i ^ 

i,3>0 

■ E Pi'(-dt)Pj'(d 8 )TN-i(t,s)\~ % ' n~ J ' 

i',j'>0 

= ]T X ~>~ m E Pi(dt)p j (-d s )r 1 (t,s)p il (-d t )p r (d s )r N - 1 (t,s). 

n,meZ i+i'=N+n-l 

Upon comparison of the coefficients of \^ n yT m for n, m > 0, we find 

T N (t,s) ]T ( L l L 2h 
0<i<N-l 

= r N (t,s) Y (4 n) (4 M ) T U 

0<i<N-l 

since = L", Lj (m) = for n, m > 0, 

= E Pi{dt)Pj(-d s )T 1 (t,s)p i '(-d t )pf{d s )T N - 1 (t,s) 

i+i' = N+n-l 
j+j'=N+m-l 

= PN+n-l{dt)PN+m-l{-ds)T\ o Tjv-1, for n, m > 0, 
leading to the statement of Theorem 6.1. □ 

7. Two-Toda symmetries and the ASV-correspondence 

Define the four vector fields 

(*i,* 2 )- = Ytf(*i,tt 2 ) := {-Nij^uNij^i) 
on the manifold of wave vectors \& = (^i, ^2) and the four vector fields 

f := XijT 

on the manifold of r-vectors. They are symmetries of the two-Toda lattice; 
i.e., they commute with the basic (t, s)-flows (see [6]). 

Theorem 7.1. There exists the Adler-Shiota-van Moerbeke- correspondence 
between symmetry vector fields on and those acting on r: 

(71) _ ( (e -„ _ 1} ^l AAe - m _ d^I). 

Proof. For iVn and -/V22, the result follows from [6] by the conversion rule 
from iVj to Nu (see footnote 13) and from (X(i, A, /ix), X(s, A, fi)) to (Xn,X2 2 ) 
(see footnote 12). For i / j, we prove for instance the following two identities: 
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For the vector fields f n = Xi2T n acting on the r- functions r n , consider the 
derivation of ^ with regard to that vector field: 

(72) (*i(*))n _ fe-^r n (t,s)\ r n (t,s) 



* n (z) ^ T n (t, S ) J e-V Tn (t, S ) 

(e' m T n (t, s))T n (t, s) - e' m T n (t, s)f n (t, s) 



= (e-*-l) 



e T 'ir„(t, s)r n (t,s) 

X 12 T n (t, s) 
T„(i,s) 



where in the above we have used the commutation of the symmetries with the 
i-flows. Considering the vector field / = X12/ (see (5.2)) acting on column 
vectors / = (f n (t, s)) ne %, we compute using (7.2) and, in the third identity, 
the relation of Corollary 3.7, with n \— > n — 1, i 1— ► t — [z^ 1 ], s 1— > s + [A -1 ], 
Z\ = z, z 2 = n, V\ = 0, V2 = A -1 , 



'V ((X 12 r(i - [Z~\ S)) n r n (t, S) - (Xi 2 T(t, 8))nTn(t - [z' 1 }, a)j 



(gj(£))n 

*in(/u) e Et i"V"T„(t, s)r„(t - [a*- 1 ], s) 



(^- 1 - Z - 1 )r n _ 1 (t-[ Z - 1 ]-[^- 1 ],s + ^- 1 ])rr l (t,s)-^ 1 r n _ 1 (t-^- 1 ],s+[\- 1 ]r n (t-[ Z - 1 ],s) 

T„(t,s)TV,(t-[/l _1 ],s) 

Tn(t, s)Tn(t - Im -1 ], S) 
= _^ Xz - ) n-l X(-8,X)X(t,z)T n -l ^ ^.^ ^ definition Qf x ^ z) 



= ~y^^2j(A by Lemma 5.4 (iii), 



*l ,.(/<) ' " ~ 
(iVl2f*l(«))„ , _ ... 

= — . . , by Proposition 6.2. 

Wl„(/U) 

thus proving 

(7.3) (flW)n = -(%*lW)„. 

Comparison of both expressions (7.2) and (7.3) for ^i(z) yields (i). 
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(n) — ^— = (Ae -1) — ■ 

Consider now the derivative of ^2 with regard to the vector field f n 
^2iT n , acting on r- functions. At first 

(7 4) ^ 2 ^ n - f e ~ V2r n+i(t,s)\ T n (t,s) 



^2n(z) \ T n (t, S ) J e-^T n+1 (t, S ) 

(e _r?2 f n+ i(i, s))r n (t, s) - (e~ m T n+1 (t, s))t n (t, s) 



= (e-^A-l)^ 



e r VT n+1 (t,s)T n (t,s) 
X 2 it x 



Acting with the vector field / = X21/ on ^2(2), and using (7.4), we apply, 
in the third identity, Corollary 3.6 with n 1— > n+1, with t s 1— > s—[v] 

and subsequently with z h A, 11 h fi^ 1 , v ^ z: 



fi , e£>H* ' [ (X 21 T(t,s-[z])) n+1 T n (t,s)-(X 21 T(.t, S )) n T n+1 (t,s-[z]) 
(^2(z)) n = V 



n+1 



c Sa i z- i -St i A i ( Z\ n ^\T n +l(t + [A" 1 ],^ - [z])T n + l(t,S - [(J,' 1 ]) 



(i\T n (t, s)r n +i(t, s - [n *]) 
1 ^ n X(-t,\)X{s,z- 1 )r n+1 (t,s) 



"(!) 

^ 7 T„(t, S ) 

= , 3>i J (A)\E'2j(^), from Lemma 5.4 (iv), 

j>n 

{N 2 lu^2{z)) n , _ ... 

= ^ — ? tt—, by Proposition 6.2, 

and thus 

(*(*))„ = (N 21u * 2 (z)) n . 

Comparison with expression (7.4) yields (ii). The proofs of the other identities 
contained in (7.1) can be done in the same style. Note that even the identities 
involving Nu and N22, which were established in [6], can be shown in this 
fashion. This ends the proof of Theorem 7.1. □ 
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8. Fredholm determinants of Christoffel-Darboux kernels 

The following theorem involves determinants of the kernels: 
K ll7n (y, z) := J2 ^u(z)^ u (y), K 2 i, n (y, z) := ^ ^* u (z)^ 2 e(y~ 1 ), 

Kn l>n 

K 12 , n (y, z) := VUz-^uiy), K 22 , n (y, z) := ]T ^(O^fo -1 ), 

Kn i>n 

already mentioned in Theorem 5.1: 

Theorem 8.1. The following holds: 

l (det (K lhn ( yi , 2j))l<iJ<fc) neZ ( det (#21,n(j/i, Z j))l<i,j<k) n 

( 8 -l) / \ / \ 

\ ( det ( R 12,n {Vi > Z j))l<i,j<k) n& ( det (^22,n (Vi > *J ) ) l<i J<fc J „, 

_ i / n*=iXn(w,«/)T rU=iX2i(y/,2/)r \ 

^ \ n*=iXi 2 (w,2/)r nJ=iX22(w,2<)r / ' 

Proof. We work out the result for the kernels if 12 and Kn. Indeed, since, 
using Lemma 5.4, 



K 12 , n (zj,yi) = ^^(Vi^uiz. 



3 j 



Kn 



V ^ r n (t,s) 
we have, using the second relation (4.1) of Theorem 4.2, 
(8.2) 

det f^ttifor 1 )*!^)) 

V <rl / l<i,j<fc 



i<ij<fc 



= (fl^ A( y - 1 )A( Z - 1 )ge E -( t ««f--f) 

T„-*(t-Ej[«r 1 ].* + EJl»r 1 ]) 



T„(t,s) 



/ fe 

* 1 



— I J^Xi 2 (zi,yi) r 
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using the computation below. Observe indeed that compounding two operators 

X12 (zi,yi) 

(Xi 2 (^i,yi)T) n = (A _1 x(yi)x(2i)^(i,^i)X(-s,yi)r) n 

= yr^r^^-^rn-iit - [zi 1 ], s + [yr 1 ]) 

yields 

(Xi 2 (22,y2)Xi2(^i,yi)r) n 

= (X 12 (z 2 ,y 2 )((y 1 z 1 r- 1 e ^(^- s ^)r n _ 1 (t-[zr 1 ],s + [yr 1 ]))n«) 

V / n 

2 



= (y^rwr 2 n e^^-^Hi - -)(i - ^) 

J=i z 2 y2 

= (n^l n^ 1 ^^^' n or 1 -*,- 1 ) n (vr 1 -^- 1 ) 

V 1 / j=l l<«<i<2 l<«<j<2 

and so on; this establishes the equality (=) in (8.2) and the ifi 2 -identity in 
(8.1). 

Since, by Lemma 5.4, 



Kn 



n 



e _ /J. _\ 



y*y i-§ T n (t, s) 

we have, using the first relation (4.1) of Theorem 4.2, 

det ]T ^ u (y t )^ u {< 

V< n ) i<i,j<k 

" ( } n 1 <„<n(y,-^) ( 1} MUrV 

7=i ' ' r n (t,a) 
1 fe 

= - TTXii(zi, yj)r. 
r i 
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In the last equality, we used the composition of the vertex operator 
X 11 (z,y)=Xl(t,y)X 1 (t,z)- 1 ^ W Sl) 



U n 1 / z — y 

several times, to yield 

k i / 1 i \ \ / \ \ / \ / ^ \ / k 



i n«(* - w) v/=i / Vo-i y ) 

(k k \ 

*+E^VE^VJ, 

thus establishing the result, for the K\\ and K\ 2 components of (8.1); the 
remaining cases are more of the same. □ 

It is also interesting to compute the Fredholm determinant of the kernel 
K = K 12 , n , namely 
(8.3) 

00 (— A) fe f f k 
det(I-XK) := 1 + E k \ J - J Ek det ( K ( x i^yj))i<i,j<kYl(p(xi,yi)dxidyi) 

over a set of the form E = E\ x E 2 C M? with regard to the measure 
p(x,y)dx dy. 

COROLLARY 8.2. The vector of Fredholm determinants (in the sense 
above) equals 

(8.4) det(/ - \K E ) = l e -^n E dxdyp(x,y)X(x,y) T 

T 

for the kernel K E = Ki 2 , n (y, z)Ie(z), with X 12 (x,y) being the corresponding 
vertex operator, given before Theorem 5.1. 

Proof. Putting the corresponding determinant obtained in Theorem 8.1 in 
the Fredholm formula (8.1), we find for a subset of the form E = Ei x E 2 C M? , 

(det(J - \K E )) neZ 

°° (— A) fc f f k 
= 1 + E fe] J ■■■ J Ek det(K(xi,yj))i<ij< k Y[{p{xi,y,i)dxidyi) 

-\) k f r i / k ^ k 



E k \ J — J Ek \ ^f[ x ( x ^yi)j T\\[p{xi,yi)dxidyi) 
E^i (~ A JJ E X(x,y)p(x,y)dxdyJ r 

}_ e -^JJ E dxdyp(x,y)X(x,y) r ^ 





1^ 1 

T 
k=0 
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9. Differential equations for vertex operators and 
a Virasoro algebra of central charge c = — 2 

Consider the vector of integrals 
(9.1) V E ■= J dx dy p(x, y)X 12 (x, y) 

over the subset E := [a, b] x [c, d] C M 2 , of the vertex operator Xi2(x, y), defined 
in (5.2), integrated over the weight 

p(x,y)dxdy := e Vl2( - x ' y) dxdy := < >: - ' ' " r/.rr///. 

Also consider the vector of operators 

with 

(9.3) j£> = (4 2 l) n& = \ (jf + (2n + k + 1) jW + n(n + l)4 0) ) neZ ; 

then the following theorem holds: 

Theorem 9.1. For all k > -1 and n > 1, 

[v fc ,(u E r] = o, 

(2) 

with the vector S k forming a Virasoro algebra of central charge c = —2: 

;4 2) , jf ] = (* - + (-2) • 

Before proving Theorem 9.1, we first need a few lemmas. For the sake of 
later investigations on matrix integrals (symmetric and symplectic), we intro- 
duce operators depending on a real parameter a > 0. So, define Heisenberg 
and Virasoro operators, depending on a, 

\ A k >° 

4V)=: k<0 and jf (a) ^j 1 *: > 

[ fc = 

together with "vector operators" acting on vectors of functions / = 
(/ra(*l,*2, • • -))n6Z, 

ifw - (4>>)„ ez 

= (j<"( Q ) + n4 0) )„ 62! md 4°' = »4 0) = ~>o*, 
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4 2) («) = (4>)) n6 



a 



E : ^(^V) = + (l " f ) ((* + 1)J? V) " W<°>) 



(f E ; : ^WjV) : + (na + (fc + 1)(1 - f )) ^(a) 

n((ra - l)a + 2) ( ) 



i+j=k 



- r 

2 fe 



(2) (2) 

Note Jj, (a) coincides for a = 1 with ' of (9.3). Given the vertex operator, 



containing a parameter a as well 

X a {u) = e^' tiU 'e 



and the "vector vertex operator," remembering \{ z ) = (■■•) z 1 , 1> z , z 2 , ■ ■ ^ 



^a(^) = A e 1 tlU e 



we prove: 



Lemma 9.2. The vertex operator X a := X a {u) and := jjf^iot) 
satisfy the relations: 



u k X a — [J^ , X a ] + X a , 
^-u k+1 (X a u n ) = £L 1 X n u n -X n u n j[ 2 l 



L « u °k,n- 



Proof. Setting X := X a {u) and using 
au-Px = \pt/3,X] and u^X 



_d_ 



X 



one immediately checks the first relation; although the following relation holds 
for all k £ Z, for brevity one checks it for k > 1, 



fc+i 



9 



<9u 



i>l 



i>l 
i + fc^O 



- E [(i-fc)ti-*,JC]«r + a E 



i>l 
-i+fe<0 



+ E lt i 

i>max(l,— fc+1) 



dt t 
d 



i>i 

-i+*>0 



5 



d_ 

dti 



dt i+ k ' 



v d 
dt k 
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E 

i>i 



it. 



d 
dt i+ k 



a \ - 

+ o E 



i+j = k 



d 2 

dtidtj 



a 
2 



E 



at,- 



+ 9 E 



X 



dtj 



+ aX 



d_ 



« fJ (2) 



+ ? E 

i 1 \x]+a^rW 
where in the last equality, one has used the identity 



+ a X J K k 



(i) 



E 

i+j=k 
»,J>1 



^ n^' W = (k — l)u k X = (k - 1) 



dt k ' . 



Finally, using the above and the first commutation relation of Lemma 9.2, one 
computes on the one hand, 



^-u k+1 (X{u)u c 
ou 



(k + l)u k X u an + an u k X u an 



(i) 



(na + k + l)u na+k X + | [jf \ Xu r ' 



- -(fc - l)u na (u k X - XJ^) + aXu na J { k 1] , 



and, on the other hand, 



jfl^Xu^-Xu^'.J, 



(2) _ 
k,n 



a 



+ [na + (k + 1) - + l)j - XJ^; 

+ 1/° + Xu na 4 0) (an + 1). 

Using again the first commutation relation of Lemma 9.2, we see that the two 
expressions are easily seen to coincide. □ 
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Proposition 9.3. Given the vector vertex operator, 

define the generators jj^ (a) by means of 
z k X a (t,z)=\^\a),X a (t,z) 



^-z k+1 X a (t,z)=[^\a),X a (t,z) 



(2) 

The S k (a)'s form a Virasoro algebra 



4 2) (a), 4 2) (a)l = (k — e)^l(a) + c [ ) <5 fe , 



12 



w/jf/i central charge 



c= 1-6 



n 



1/2 



a 



-1/2N 



TTie J^(a) /orm a Heisenberg algebra, interacting with Virasoro, as follows: 



4 2) (a),J«(a) 



A; 



Q: 



4, 



1 1 



a 2 



Proof. The proof follows from Lemma 9.2 and an explicit computation for 
the central charge. □ 

Proposition 9.4. The vertex operator, defined in (5.2), 

X 12 (u,v) = K- 1 X{-s,v)X{t,u) X {u)x{v) 

leads to a Virasoro algebra of central charge c = —2, 

—u k+1 X 12 (u,v) = [4 2) ,Xi 2 («,t;)] , 



with generators (in t) 
4 2) :=4 2) (a) 



a=l 



J, 



(2) 



k,n 



raGZ 



i(jf ) + (2n + fe + l)4 1) +n(n + l)jf) 



Similarly, the involution u <-> -u, f —s leads to the same Virasoro algebra 
in s, with same central charge. 
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Proof. Noticing that one piece of Xi 2 (t, s;u,v) is precisely X Q (t;«) for 
a = 1, we apply Proposition 9.3 for a = 1: 



<9u 



u fc+1 X 12 (i,s;u,t;) 



d_ 

du 

'4 2) (l),Xi (*;«)] X(-s,v) X (v) 
ir( 2 ) 



= [4 ; ^ 12 ]. 

The central charge c = — 2 is obtained by setting a = 1 in the general formula 
for c in Proposition 9.3. □ 

Proof of Theorem 9.1. We consider the vector of operators given by the 
double integral of a vertex operator, 

(9.4) j\ x jy y JL ( x ^X 12 (x,y)p(x,y)) , 

thought of as acting on a column of functions F(t, s, c). On the one hand, the 
integral (9.4) equals 



(9.5) 



= x k+l 



,Vi2(x,y) 



J dyX 12 (x,y)e 

bk+1 §b + ak+1 ~L) Ia dx Ic 1 d y^ x iy) eVl2{x,v) ■ 



x=a 
d 



On the other hand, by Proposition 9.4, it equals 

rd 



(9.6) 



+ 



dy J dx (j^x k+1 X 12 (x,y)^j p{x,y) 

( dy ( dxX 12 (x,y)x k+1 ( icijX l ~ V \ p(x,y) 
Jc Ja \ij>i J 



id 



dy / dxX 12 (x,y)e Vl2{x > y) 



ij>l dci+k ' j 

where the derivations d/da, d/db, d/dc act on the operator only and not on the 
function F(t,s,c); note that for an operator A, {dA)F = [d,A\F. Comparing 
the two ways (9.5) and (9.6) of computing (9.4), we obtain [Vfc,Ue] = 0, with 
Vfc as defined in the beginning of this section. The rest follows from the next 
argument: 

n 

[V,,U£] = ^U^[V f ^Ul" 1 = 0. 
k=i 

Incidentally, in view of Corollary 8.2, it also implies that \Vg, e _AUi5 ] = 0, since 



-AU F 



Eoo 
n=0 



(-xy 



it, 



i! E' 



□ 
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10. Vertex representation of probabilities 
and Virasoro constraints 

Consider a weight p(y, z)dy dz := pt yS (y, z ) '■= e Vt ' a ^ y ' z ^dy dz on M. 2 , with 
po = e v °, where 

00 oo 

(10.1) V t , s {y, z) := cyz + £ Uv* ~ s ^ ■ 

1 l 

Given the space of Hermitean matrices TCn, and given 

spectrum M\ = {x±, . . . , xjy} and 

spectrum M 2 = {yi, ■ ■ ■ , y./v}, with Mi, M 2 G Hat, 

we define, for a set E C M 2 , 

Wjv,js = {(Mi, M 2 ) G ?4 with all (x k , y t ) G £}. 

Consider the product Haar measure dM\dM 2 on the product space H 2 N , with 
each dMi, decomposed into its radial part and its angular part, as in (0.16). 
Also define the probability measure 



dM 1 dM 2 e TlVt ^ Ml ' M ^ 
JJ H 2 dM 1 dM 2 e T " v ^ M ^ M ^' 



(2) 

Recall from (9.3), the definition of the vector , also define another one 

(10.2) = (4 2 ) )neZ = I(jf) + (2n + fc + l)J«+n(n+l)4 0) ) ngZ , 

ft 2) = (j£)n eZ = ^(jf ) + (2n + fc + l)4 1) +n(n + l)4 0) ) neZ . 
Given the disjoint union 

(10.3) E = E 1 xE 2 := U\ =x [a 2i -i , a 2i \ x Uf =1 [6 2i _i, 6 2i ] C M 2 , 
define the following integral: 

(10.4) U E -= Xi 2 (x,y)p (x,y)dxdy, 



IE 

of the vertex operator X12, defined in (5.2). The main theorem of this section 
is: 



14 J* = 



1,2. 
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Theorem 10.1. Given the set E, as in (10.3), the probability 
P(all M\- eigenvalues € E\ and all M 2 - eigenvalues € E2) 

JJ H 2 dM 1 dM 2 e TrV ^ M ^ T E 
1 ' H Hl dM l dM 2 e^A M uM 2 ) ■ Tn 

is a ratio of two r -functions and r n , such that 

rf = ((U £ )V) n . 

Moreover, r n and r,f satisfy the partial differential equations, labeled for 
k > -I, 

(10.6) ^-jy^^ + tf^ = 

Remark. Whenever some a% or foj = 00, we must interpret: or 
= 0; in particular r n satisfies the same equations, but without the 
boundary terms. 

The following proposition is due to [12], [8], [9]: 

Proposition 10.2. 



(10.7) / dU<f*«>*> T = (2?r) , 2 ^T'^ . 

JUM n\ A(x)A(y) 



I 2 _ 

lu{n) ~~ ~ n\ " A(x)A(y) 

Proposition 10.3. For E = E 1 x E 2 C M 2 , i/ie following holds: 
(10.8) / / e cT ^ M ^e T ^7(^i-s^) dMidM2 



r r N 

/ / ]\(dx k dy k eV£ite*l-«vV-**M)A N (x)A N (y). 
J JeN k=i 



Proof. Consider a symmetric function f(x, y) := f(x\, . . . , xjv; yi, ■ ■ ■ , I/at) 
in 2/1, . . . , yjy for given xi, . . . , xn] then we have, using the skew-symmetry of 
the Vandermonde Ajv(y), 

/ / A 7V (x)Aj V (y)/(x,y)det(e :r *)i<i,<iv^y 

A n (x) V (-irA w (y)/(x,y)e S ^«^dy 
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&n(x) (- 1 ) CTA A r (y CT -i(i),---,y CT -i(iv)) 

aen N 

• f(x; Va-Hi), ■ ■ ■ , Va-^N))^ 1 X ^dxdy 



R 2N 



^n(x) ]T (" i r(- i r ) A N ( yi ,...,y N )f(x,y)e^ x ^dxdy 



N 



N\ f [ A N (x)A N (y)f(x,y)e^ x ^dxdy. 

J JR 2N 



The function 



N N 

f(x,y) =Y[iE 1 xE 2 (xi,yi) = Y[l El ( x i) I E 2 (yi) 
i i 

has the desired symmetry property, so that, using (10.7), one computes 
f f e cT ^ M ^e Tr ^^ M l-^) d M 1 dM 2 

■* ^N,E 

N N N 

= 11 n^II^A^(:r)A^(y) f[ e^'^-^) 

J JeN i l k=i 

■ [ [ d\J x dV 2 e c ^ v ^ v ^ 

J JU(N) 2 

r r N N N 

= 11 n^n^w^^n^ 1 ^- 3 ^ 

j jen i i k=i 

■ 1 1 dU ld U 2 e cTrxU ^, 

J JU(N) 2 



U(N) 2 

substituting U 2 for UJU2, 



N N 



c N I I II dxt 1] dy i A 2 N (x)A 2 N (y)f{x, y) 



tA A N (x)A N (y) 



= c 'n II \{d Xl \{dy t A N (x)A N (y) f[ e ^=^- s ^ +c 
J JeN i l k=i 

where the last identity follows from the previous calculation. □ 
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Consider now a more general weight p(y, z)dydz := pt, s (y, z)dydz :- 
e Vt - s ( y > z ^dydz on R 2 , with p = e v °, where 



CO oo oo oo 



(10.9) V tjS {y, z) := V (y, z)+J2 Uv*-J2 ^ = E C ^ j +E s 

1 1 i,j>l 1 1 

with arbitrary Vo and the inner product with regard to a subset Eci 2 

(10.10) (f,g)E= I dydzp t , s (y,z)f(y)g(z). 

Je 

Given the moment matrix (over E), 

(10.11) m n (t,s,c) =: (pij) <ij<n-i = ({y l ,z j ) E ) <ij<n-i, 
according to [2], [3], the Borel decomposition of the semi-infinite matrixE! 

nioo = S± S 2 with Si 6 £>-oo,o, S 2 G V 0)OO 

with Si having l's on the diagonal, and S 2 having hfs on the diagonal, leads to 
two strings (p^\y) , p^ 2 \z)) of monic polynomials in one variable (dependent 
on E), constructed, in terms of the character x(z) = ( zn )n&,n>o, as follows: 

(10.12) p«(y) =: Six(y), P {2 \z) =■ h{S^) T x{z). 

We call these two sequences hi- orthogonal polynomials; in fact, according to [3] 
the Borel decomposition of = S± 1 S 2 above is equivalent to the "orthogo- 
nality" relations of the polynomials: 

(10.13) {p { n\p i u})E = 5 n>m h n . 

The matrices 

Li := SiASi 1 , and L 2 := S 2 A J S? 1 , 
interact with the vector of string orthogonal polynomials, as follows: 

(10.14) Lip^Xy) = yp {1 \y), hLj h' 1 p {2) (z) = zp^{z). 
Also define vectors ^>i and as follows: 

(10.15) #i(z) := e^ zk p^\z) and %(z) := e^^'Vy 2 ^" 1 ) 

= e^ zk Six(z) = e-*">*- h {SZ 1 ) T x{z- 1 ). 



5 I>fc ^ (k < £ £ Z) denotes the set of band matrices with zeros outside the strip (k, 
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As a function of (t,s), the couple L := [L\,L-2) satisfies the two-Toda lattice 
equations (2.3), and and ^ satisfy the equations (2.7); remember that L, 
and *2 au depend on E. 
Moreover, according to [2], Theorem 3.4, the determinant of the moment 
matrix can be expressed as a 2n-uple integral over E n C R 2n : 

(10-16) 

n\detm n (t,s,c) =11 A n (u)A n (v) TT (e Vt > s( - Uk ' v ^du k dv k ) 

= n\det(E n (t) moD (0,0,c)E n (-s) T 



n-1 



= ]Jhi(t,s,c) 





where E n (t) := (the first n rows of e^ 1 * nA ™) is a matrix of Schur polynomials 
p n {t). Also T n (t,s,c) is a r-function with regard to t and s. Note that for 
Vq = cxy and E = E\ x _E 2 , this integral is precisely the one obtained in 
Proposition 10.3. 

Proposition 10.4. Given the bi- orthogonal polynomials {p^\ pj^) f or 
a general weight po = e v ° on M 2 , the kernel defined in terms of (10.15), 

(10.17) 

K(y,y';z,z'):=K n (y,z') := £ ^ 1;k (y)^ k (z - 1 ) 

0<k<n 

= E e^p^\y)h k ^\ Z ')e-^r^\ 

0<k<n 

defines a projector; i.e., it has the reproducing property with regard to the 
measure podz dz': 

I \ K(y, y'; z, z')K(z, z'\ u, u')p (z, z')dz dz' = K(y, y'; u, v!) 

and 

(10.18) // K(z,z';z,z')p {z,z')dzdz' = n. 



Proof. Using the explicit expression (10.17), using the fact that pt <s (y,z) 
= po(y, z)e^' 1 (^-s^ 1 ) i n the second equality, and using the orthogonality 
relation (10.13), in the third equality, one computes 
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j j K(y,y'; z, z')K(z, z';u,u')po(z, z')dzdz' 

J JR 2 

lii E ^\y)K^\z') 

J JR \0<k<n 

]T <^ zi p£\z)hJ 1 p? ) (u')e-* 9ivfi ] Po (z,z')dzdz' 

)<i<n 

= // E ^Vi%)/^fV) 



0<fc/<n 

-S^ (2), ,x (1), _l 



•e-^ s ^ p%> (z') P y> {z)h?fH,,{z, z')dz dz> 



0<k<n 

and 



/ / K(z, z' ; z, z')po(z, z')dz dz' 

J JR 2 

= // i^n^, z')p (z,z')dzdz' 

J JR2 



/ / 

0<fc<n' 

n. □ 



= E I / R ^\z)p[ 2) (z')h^p t Az,z')dzdz' 



Consider, for E G M 2 , the vectors (see (10.16)) 
(10.19) r E ={r r f) := Iff f[(dx k dy k p t<s (x k ,y k ))A n (x)A n (y) 

r ■= K) n > :=(r f f) 



" k=l / n>0 



'n>0 



with po as hi (10.9), and the vector of operators U^, defined in (10.4), but for 
the weight po as in (10.9). Given the set E C M 2 , define, in accordance with 
(9.3): 



(10.20) V, := -E^ +1 ^+i 2) +E-^ 



(10.21) 



v 4 , + 



i=l i,j>l 

We now state: 
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Proposition 10.5. For E = E\ x E 2 C M 2 ,@ we have 
(10.22) rf = {{V E ) n r) n . 

Proof. In what follows, we use the monic bi-orthogonal polynomials p^' , pp , 
defined by pt, s ( x ,y) on M 2 ; therefore the hi(t,s,c) are the M? inner products. 
We first compute, using (10.16) for E = M 2 and Proposition 10.3, and remem- 
bering notation (10.4), and formulae (10.16) and (10.17): 

E ( n ~l \ p p n 

T~ = II^ 1 // Y\_( dx kdy k pt,s(x k ,y k ))A n (x)A. n (y) 
\ o / *' J En fe=l 

(Ti—i \ r r n 
1 J J J IX^ 35 * 5 ^' 9 *' 8 ^*' yfe )) ^(Pi-iC^j))!^*^'^" det (Pi-i(yj))i<i.j<» 

p p n / n 

= // Y\_{dx k dy k p (x k ,y k ))det e SM! *> ^p^^/CVPi-ife 
J «" k=i \ i=i 

Y\_{dx k dy k p (x k ,y k ))Aet\ ^ ^liOfc)*^^ 1 ) I 

ft = l \0<i<n-l / 

/ / J\(Podx k dy k ) det (Knixi^ye))^ e < n 
J Je " k=i 

n ( \ 71 \ 

Y\{podx k dy k ) - ]^[Xi2(a; fc ,y fc )r J , using (8.1), 
\ fc=i 

X 12 (x,y)p (x,y)dxdy 



establishing (10.22). □ 

Proof of Theorem 10.1. In [2] (see for instance the introduction), we have 
shown that the vector r = t k2 , which is independent of the a^'s and frj's, 
satisfies the infinite set of equations, for k > — 1, 

Vfcr = (4 2) + £ icy-J?_ ) T = 



VfcT = (j1 2) + 2 j Cij ._^_)r = 



16 ((Ufi)™i")n means the n th component of the vector (Ub)°t. 
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According to Theorem 9.1 , we also have \V k , (U^;)™] = and thus 

= [V fe , (V E ) n ]r = Y k (V E ) n r - (V E ) n Y k T = Y k (V E ) n r; 

taking the n th component, we find (¥ k (U E ) n T) n = and similarly with Y k 
replaced by V k . Since ((U E ) n T) n = t£ by (10.22), this leads to: 

(^) {-i^^ + J^E^T^-}^ = 




i=i ij>i 

But, by p. 285 of [2] and by Theorem 6.1, 
r)T E n ~ 1 

(10.24) = T*J2( L i L $h = Pa+ n-i(dt)pp + n-i(-d s )rF o rti- 

Remember, one is really interested in the probability, expressed by r-functions 
(see (10.8) and (10.16)), 

JJ H 2 dM 1 dM 2 e TrVt ^ T E 

for po = e Vo = e cxy ; thus, we must set all Cij = 0, but c = en; this leads to the 
statements (10.6), ending the proof of Theorem 10.1. □ 

11. PDE's for the joint statistics of the spectra 
of Gaussian coupled random matrices 

Consider the Gaussian probability measure 

(11.1) c n dM 1 dM 2 e-^ M ? +M '- 2cM ' M2 \ 

defined over the space of Hermitean matricex = H n x H n , with a coupling 
constant c. Consider the joint probability 

(11.2) P n (E) := P(all (Mi-eigenvalues) G E 1: (M 2 -eigenvalues)e E 2 ) 

for a set of the form E = E\ x E 2 := U[ =1 [a2i-i, a 2 i] x Uf =1 [&2i-i, &2i] C R 2 . 
Before stating the theorem, we remind the reader of the differential operators 
A k , B k , depending on the boundary points of E and the coupling constant c: 
(11.3) 

1 fr q s q X ! ( ' d ' 9 \ 
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1 + c 2 , r , „ , 2c 



they form a Lie algebra parametrized by c: 

(11.4) [A 1 ,B 1 ] = [A 1 ,A 2 ] = ^^A 1 [A 2 ,B 1 ] = -^A 1 

1 — c z 1 — cr 

1 — c z 1 — c z 

We now prove Theorem 0.5, as announced in the introduction: 

Theorem 11.1 (Gaussian probability). The joint statistics (11.2) satis- 
fies the nonlinear third- order partial differential equatiow-H (F n := — log P n (E)): 
(11.5) 

B 2 A x F n , B x AxF n + -^—) - {A 2 B x F n , A x B x F n + -^—] = 0. 



Remark. When E\ = E 2 , equation (11.5) is trivially satisfied. 
Proof. From (10.5), it clearly follows that 

T n(t> s > c ij) 



Pn(E) 



T n s i c ij) £ 

where t® is an integral over E n C M. 2n , i.e., (x, y) E E™ x E% = -E n , 
(11.6) 

Tn^^^ij) = // dxdyA n (x)A n (y) 

J JE n 

jl -|(4+2/I- 2ca; fc»)+E l ^i(**4- s i4)+E i,j>a 

fc=l 

and where C denotes the locus 

C = {U = Si = 0, en = c and all other dj = 0} . 
Observe the following involution on r^(t, s) = T n (t, s, a, b, c): 

(11.7) r n (-s, -t,b,a,c) = T n (t,s,a,b,c), 

implying for the Ai, Bi, defined in (11.3) and Vi, Wj, defined below: 

Ai ~ (-1)%, V t <-» (-1) ,: W« . 

In view of (11.6), we write down the Virasoro equations (10.23) for t^, 
but with the shifts t 2 i— ► — 1 + 1 2 , s 2 i— ► 1 + s 2 - It is convenient to consider new 
Virasoro generators and Wk, such that 

(11.8) V ^ = ^ and W ^ = ± ^ 



17 in terms of the Wronskian {/, g}x = Xf.g — f.Xg, with regard to a first-order differential 
operator X. 



970 



M. ADLER AND P. VAN MOERBEKE 



namely, in terms of (10.20) and (10.21): 
(11.9) 

l 



Vi = 



, fV-i + cV.i) 

Cr — 1 

d n(ti — csi) 
Tin + c 2 - 1 



/ 



+ 



c 2 - 1 
1 



E., d d . v— v .. d d . 
i(ti-KL 1" c «i » )+ > <Wa 1" J c a ) 



\ 



i,3>l 



/ 



Wi = — L.-( c V-i+V-i) 
1 — c 2 

9 n(cti — si) 
dsi c 2 — 1 



c 2 - 1 



9si_i ^ da-i,3 ddj-i' 

i,j>l 



V 2 := Vo-c 



9c 



N _a_ n(n + 1) 



E 



»>i 



M>1 



W 2 := Vo-c— 



<9s2 ^-^ * A.<?„ 



9 n(n+l) 



9s, 



E ■' r;j 7h~ 



i,j>l 



With this new notation, and by virtue of (10.23) and (10.2), the r ra 's satisfy 
for all n > 1: 

(11.10) A k T n = VkT n and S fc r n = W fc r n , fc = l,2. 

In particular, on the locus £, we have from (11.9), 

B\T n 



A\T n 
A 2 T n 



dti 



dsi 



( d_ , n(n+l) \ 

\ dt 2 + 2 ; 



and so 
(11.11) 



Jr lo s r ™ 



Ai log r n 



»(»+i) \ 

V<9s 2 + 2 J r ™ 



2 9s 2 10 & T ™ 



#i logr n 



#2 logr r , 



n(n+l) 
2 
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Using [Si,Vi] 



0, we have 



BiAlTr, 



V x B x r n 
d f d 

( d 2 



ct\ — S\ 



dt\ds] 



+ 



nc 



1-c 2 



and so, on the locus C, 

d 2 



(11.12) 



Using [B 2 ,Vi] 



dt\ds 



logr n 



B\A\ logr n + 



nc 



0, we see that 



B 2 A\T n 



B 2 V X T n 

c 

ViB 2 T n \ c 
ViW 2 T n \ c 

d f d n(n + 1) 
dh \ds2~ + 2 



d 2 n(n + 1) d 



dt\ds 2 



dti 



and so, on £, we haveE 
(11.13) 



d 2 



dt\ds' A 



■ log T r , 



B 2 A\ logr n . 



Setting (11.11), (11.12), (11.13) into the formula of Proposition 3.3 (for k = 2, 
as spelled out in Lemma 3.4) and its dual, namely 

9 .iog rn+1 - a5fe logTn 



(11.14) 



3 \og Tn+l ~ ^fc l0gr ' 



ds i r, 



and 



n-l 



dt\ds 



Using the following relation for non-commutative operators X and Y 
XY\ogf=j 7 (fXYf-Xf Yf). 
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one is led to an expression for B\ log and, using the involution (11.7), a 
dual expression for A\ log : 

T n +1 A 2 B\ log T n 



(11.15) -Ai log 

-Hi log 



r n _i ^li^i log r n + 

T n+ l B 2 Al log T n 



r n _i Si^lilogr n + 

Upon taking A\ of the second expression, subtracting from it B\ of the 
first one and using [.Ai,£>i] = 0, one finds the following identity 

B 2 A\ log T n _ ^2^1 log T w = 

1 SUilogr n + ^ 1 ^ 1 l gr n + ^"' 

This difference amounts to the equality of two Wronskians (G n := ^ logr n ): 
(11.16) 

B 2 A 1 G n , B\A\G n + 1 = (.A 2 £iG n , A x B x G n + 



c 1 J Ai L c 1 J Si 

Because of the fact that 

logP n (£) = log(r n (E)/r n (M 2 )) = logr n (£) - logr n (M 2 ), 

together with the fact that ^4ir n (IR 2 ) = BiT n (M. 2 ) = 0, we have that 
F n (E) := \ log P n {E) satisfies (11.16) as well, thus leading to (11.5). □ 

Remark. For small n, the equation (11.5) for r n = detm n can be checked, 
using the explicit moment matrix m n = (fJ,ij)o<i,j<n-i, where 

= f dx f dyx l y ] e-^ x2+y2 - 2cxy \ 

J Ei JE 2 

It suffices to compute the action of Ai and B; L on fiij, namely 
(11.17) 

ifii-ij + cjfiij-i 

A\Hij = H ' 2 _ | — '■ A 2 fiij = {i+ 1)/%- — Hi+2,j 

Bifiij = -tMj+i + i_ c 2 B 2M»j = U + !)/%■ - Mi,i+2 , 

and check equation (11.15), at least for small n. 



12. Coupled random matrices with 
the Laguerre statistics 

Consider the Laguerre probability measure 
(12.1) c n dM 1 dM 2 e Tl ^ Ml+a log Ml - M2+a log M ^MiM 2 )^ 
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defined over the space of Hermitean matricex TC^ = 7i n x TL n , with a coupling 
constant c. Consider the joint probability 

(12.2) P n {E) := P(all (Mi -eigenvalues) G E u (M 2 -eigenvalues)e E 2 ) 

for a set of the form E = E\ x E 2 := U^ =1 [a 2 i-i , a 2 i] x U* =1 [&2i-i > &2i] C M 2 . 
Before stating the theorem, we remind the reader of the differential operators 
Ak, Bk, depending on the boundary points of E and the coupling constant c: 

(12-3) 

V s 

-4 2 = y^os?-a— + (n+l + a)Ai -c- , B 2 = y^b 2 j — + (n + l+a)Bi - — . 

Note that A 2 and #2 acting on r n depends on the index n. 

Theorem 12.1 (Laguerre distibution). The joint statistics (11.2), namely 
P n {E) = T n (E)/T n (M + ), is a ratio of two functions, each satisfying the non- 
linear third-order partial differential equation^ (G n := log r„(i?)): 

(12.4) 

{(B2A1 + cJ-)G n , BxAxGrX - {{A 2 B X + c^-)G n , A^gA = . 

I <9ci2 J Ai I <9c 2 i J Si 

Remark. Equation (12.4) is actually an inductive set of equations with re- 
gard to n, since it contains derivatives of the form dr n {E) / dc 2 \ and dr n {E) / dc\ 2 . 
The point is that, according to (10.24), these derivatives can be expressed in 
terms of (t, s)-derivatives of the expression T n -\(E) in (12.5) below; to be 
precise, 

8r n (E) 



Do 



21 



Pn+l(d t )pn(-d s )Ti(E) O T n -i(E) 



t=s=Cij=0,cn=c 



E fif jPi >(-d t ) Pjl (d s )T n ^(E) 

j+j'=n 
i,i',3,f>0 



t=s=Cij=0,cn=c 



The t, s-partials of T n -\(E) can then be expressed again in terms the operators 
Ai,Bi applied to r n _i, etc... . 

This result hinges on knowing, as before, the Virasoro constraints for the 
(t, s) deformations of the matrix integral (12.1). Unlike the Gaussian case, 



19 in terms of the Wronskian {/, g}x = Xf.g — f.Xg, with regard to a first-order differential 
operator X. 
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which could be obtained by merely shifting the time, we invoke here a method, 
due to [5], of representing the matrix integral by vertex operators acting on a 
vacuum vector. 

Proposition 12.2. The integral 
(12 ' 5) 

r n {E) = If dxdyA n (x)A n (y) f[ p(x k )p(y k )e^ ? = l{ti ^ s ^) + ^>i c -^, 
J jEn fc =i 
with the weights p and p satisfying 

(12.6) -P- = l = ^ and -9 Ei>o ft* 



P f Ei>O a i zt P f Ei>0«^ 

satisfies the following Virasoro equations for k > — 1: 
(12.7) 

(E^ +1 /(ao^--E(-(^ + E ^db'-^OK 

\ 1 i>0 v m,t>l ' / 

(e^<<-E(*«<+ E < ~adb»-* J 2«)) T - a 

\ 1 i>0 V m,£>l 7 / 

Proo/. See [5]. 

Proof of Theorem 12.1. Since 

p(z) = /5(z) = z a e~ z , 

with 

— — = and f(z) = z, 

p z M ; ' 

we have /?o = —a, /3i = 1, ao = 0, ai = 1 and all remaining «j and = 0. 
Therefore, for all k > 0, 



\ 1 OT i \ y>, 



where 1^,1^,1^,1^ are defined in (10.2) and in the formulas following 
Theorem 3.1. Guided by the same principle as (11.8), one redefines 

d d ^ d n(n + 2a + l) 
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^. d d x . 

V2 = 5>^7-^ + E 



lc ij— + (n + a + l)(Vi + 



0tr 



and similarly for Wi, using the map s <-> — i and icijd/dck+ij «-> jcijd/dcij + k- 
Here the involution acts as follows: 

A <-» Sj, V^W,. 

Then r n satisfies for = 1, 2 

(12.9) A^n = V fe r n and S fc r n = W fc r n . 

Evaluating AkT n , BkT n , A\B\T n , A 2 B\T n along the locus £, using the 
commutation relation [A2, Wi]\c = — c a^" an< ^ (12.9), and setting d n = 
n(n + 2a + l)/2, one checks 



9r n 



(^i-d n )r n | £ — 



c 



dsidti 
d 2 r n 



= -(Ai- d n ) (Si - dn) r n \ c 



ds\dt 2 



(Bi - d n )r n \ c 



d 



(A 2 - (n + a + l)d n ) {Bi - d n ) r n ~ c 'g^ T n 



and so, 



d 2 log T n 



d log T„ 



ds\dt\ 



-A\B\ logr n | £ and 



Si logr n | £ - 

<9 2 logT n 



n(n + 2a + l) 



ds\dt 2 



d 

A 2 Bi + c- — ) log r n 



Setting these expressions into equations (11.14), we find an equation and its 
dual: 

_ Alog -±i +( , B+2a+1) _ ^Y;^ )logT " 

r n _i ^i/5ilogr n 
_ Bllog Z»±l + (2 „ + 2a+1) _ ^-4. )l0gT " ■ 

Finally, upon subtracting £>i of the first from Ai of the second, one is led to: 

_ (A 2 B 1+ c^)logr n (S 2 A + c^)logr ra 
£>1 ^ — -Ai — — - A — = U, 



A\Bx logr„ 
ending the proof of Theorem 12.1 



B X A\ logr n 
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